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FRS  SCOPING  STUDY  D:  ZPIMP  PREDICTIONS  FOR  RADUTION 
YIELDS  ON  A  JUPITER-CLASS  GENERATOR 

I.  Introduction 

The  Defense  Nuclear  Agency  (DNA)  has  long  supported  plasma  radiation-source  (PRS) 
development  for  assessing  radiation  hardening  of  various  defense  and  communication  systems. 
Existing  pulsed  power  generators  sudi  as  DOUBLE  EAGLE.  BLACKJACK  S.  PHOENIX,  and 
SATURN  deliver  from  3  to  10  MA  of  current  to  an  imploding  z  pinch  and  can  produce  25  -  75  kJ 
of  A£  K-shell  radiation  (>  1.6  keVirfiotoos)[R^.l].  The  DECADE  generator  under  constructimi 
is  designed  to  deliver  ~20  MA  of  current  to  the  front  end.  With  a  Zrfwch  load,  DECADE  is 
predicted  to  produce  ~300  kJ  of  argon  K-shell  emission  (>  3.1  keV  photons)  [Ref.2].  Mtii  the 
cessation  of  underground  testing  die  next  generation  pulsed  power  driver.will  be  required  to  produce 
mega-joules  of  radiation  in  tiie  0- 15  keV  range,  and  several  hundred  Idlo-jouks  between  30  and 
60  keV.  The  JUPTIER  Design  Option  Study  Team  (JDOST)  was  formed  in  the  summer  of  1993 
to  develop  circuit  models  for  tins  next  generation  mKJiine.  R  is  anticqMted  that  tiie  JUPTIER 
madiine  wQl  provide  ~60  MA  peak  load  current  In  order  to  determine  tiie  potential  radiation 
yidds  from  sudi  a  madune,  DNA  has  duutered  tile  Plasma  Physics  Division  at  the  Nanral  Research 
Laboratory  to  undertake  a  PRS  soopmg  study.  This  is  me  of  several  pliers  reporting  on  tiie  results 
and  condusicms  of  that  study. 

In  die  precursor  pqier  [Re£3]  to  tiiis  rqxxt  tiie  Whitn^'Thornhill  scaling  law  [RefB.4  and  5] 
was  gmeralized  to  account  f(v  two  problems  diaracteristic  of  inqiloding  loads  on  siqier  hi^  current 
genaatms:  0)  lari^  initial  radii,  and  0i)  weak  inqilosions  which  batdy  tbermalize  lugh  atomic 
number  material  like  Kr  into  the  K-shell  ionimtion  stage.  The  resultant  algdmic  relations  for  the 
radiation  yield,  termed  the  J-scaling  law,  depead  on  simple  prcqierties  of  tiie  inqilosion:  atomic 
number  Z,  mass  loading  Af,  pinch  length  f,  final  kinetic  oieigy  Ek*  and  stagnation  radius  Rf. 
The  overriding  utility  of  the  scaling  qifnoach  lies  in  the  rapid  and  efficient  investigation  of  large 
ranges  in  load  and  madiine  paramder  space.  Refnmce  3  used  the  J-scaling  law  in  conjunction 
with  citcuit  models  to  estiirutte  yields  <»  tiie  two  leading  design  candidates  for  JUPITER:  the  Umar 
Inductive  voltage  Addo’  (LIA)  and  the  Inductive  Energy  Store  (lES)  designs.  The  peak  yields  for 
the  standard  point  designs  of  botii  machines  were  fcMind  to  range  from  ~7  MJ  for  Ar,  throu^  ~3 
Ml  for  Kr,  to  MJ  for  Jfo.  Furthermore,  the  J-scaling  law  was  used  to  surv^  the  Ar,  Kr, 
and  Xe  yield  dqiendmioe  on  pinch  length,  front  end  inductance,  and  switdi  behavior  for  various 
modifications  of  the  standard  designs. 

The  ease  of  triplication  of  the  J-scaling  law  belies  the  large  investment  in  its  development. 
The  precursor  to  the  J-scaling  law  was  based  on  and  tested  agtunst  Ai  experiments  on  existing 
marJiiiiftB  [Refs.4].  In  this  regime  one  finds  good  agreement  with  the  data  (Ref.l].  However,  Ref.3 
noted  serious  issues  in  ^tending  the  K-shell  scaling  law  to  JUPITER  load  configurations.  First, 
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it  is  highly  unclear  whether  the  Z  scaling,  derived  for  low  atomic  number  elements,  applies  to  the 
umization  dynamics  of  kiypton  and  xenon,  where  the  L-and  M-shell  states  can  act  as  a  major  eneigy 
sink.  Second,  the  stagflation  physics  at  ~60  MA  on  JUPITER  may  present  significant  differences 
from  the  case  at  ~10  MA,  where  the  scaling  was  tested.  Third,  the  optimized  yields  on  JUPITER 
are  found  at  large  radii  but  the  effects  of  tmiahilittff,  asymmetries,  and  turbulence  are  unknown  in 
thisrq^ime. 

The  goal  of  die  present  rqxxt  is  to  irqxove  iqxm  die  scaling  law  yield  estimates  by  addressing 
in  more  detail  die  first  two  issues:  ionization  and  stagnation  physics.  A  simple  1-D  simulation 
code,  2PIMP,  will  be  enqdoyed  to  follow  die  z  pindi  tfynamics  tm  a  JUFTIER-class  generator 
througjiin:ylosi<Mi.diermali7ation,  and  subsequent  bounce,  while  coupled  to  an  atomic  irfiysics  and 
radiatkm  tranqxxt  fflodeL  This  will  provide  an  improved  assessment  of  the  PRS  cqiabilides  on 
JUFTIER,  lead  to  a  better  understanding  of  (he  unique  stagnation  id^sics  at  siqier  hi|^  currents, 
and,  because  of  the  ciqMbility  of  qiectnl  resolution  in  a  full  radiation  calCTilation,  allow  a  direct 
copqiatison  with  stated  radiadon  requirements  in  various  energy  bins.  It  should  not  be  presumed 
thafaeding  law  ealrMlarinn«  am  jfTdewant  in  fhftfteeofmomcnmplex  modeling.  AswillbepOUlled 
out  dnon^KNit  dus  papers  scaling  law  estimates  wifi  be  used  to  eacamine  die  gross  dqiendence  of 
die  yield  in  (he  load  mass  -  initial  radius  plane.  This  allows  die  ZPIMP  studies  to  focus  on  a 
sdectedr^tm  of  parameter  space  widi  a  limited  set  of  tuns  to  determine  die  qpdmal  coodidons. 
As  the  ZPIMP  code  {xovides  miniimJ  spadal  resoludon  by  treating  the  plasma  shell  as  a  sin^ 
zone,  it  too  wiU  iMed  to  be  rqilaced  widi  more  cooqilex  muld-zone  and  VD  simuladons.  One 
should  view  the  scaling  law  followed  by  the  ZPIMP  calculadons  as  the  b^inning  in  a  sequence  of 
steps  leading  to  more  refined,  coiiq)lex,  and  more  time  consuming,  3rield  estimates. 

Section  n  contains  an  overview  of  the  conqwnents  in  the  ZITMP  (Z-Pinch  IMPlosion) 
simulation  code.  The  equations  desoibing  the  dynamics  of  the  plasma  shell  as  a  single  zmie 
are  summarized  along  with  the  magnetic  induction  equatimi,  the  lumped  circuit  model  driving  the 
pinch,  the  atomic  physics  models,  and  the  radiation  transport  technique.  A  low  mass  core  zone 
inside  the  plasma  shell  is  also  modeled.  Since  the  code  takes  a  Lagrangian  viewpoint,  an  artificial 
viscosity,  ,  is  defined  to  effectivdy  soften  the  imploaon  [Ref .5] .  In  Sectimi  ni  a  value  f m  is 
detmnined  by  matfibing  the  K-shell  yield  from  a  large  set  of  SATURN  aluminum  experiments. 
The  resultant  code  is  thm  used  to  investigate  krypton  implosions  on  a  JUPITER-class  generator  in 
Section  TV.  The  JUPTIER  model  in  the  present  paper  is  a  Thevenin  equivalent  circuit  of  an  early 
LIA  design.  For  the  9«Mappropriate  to  SATURN,  JUPITER  promises  poor  yields  in  the  higher 
photon  ermgy  bins.  To  enhance  the  yields  the  pinch  coixq)ression  at  stagnation  is  increased  by 
arbitrarily  lowering  q^.  Section  V  presents  a  list  of  radiation  fluence  test  requirements  in  various 
energy  bins  for  JUITTER  and  converts  them  to  source  yields  subject  to  uniformity  criteria  and 
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debris  standoff  distances.  Using  2PIMP  with  ttw  above  mentioned  adder  circuit,  a  search  over 
load  mass  and  initial  radius  is  performed  to  optimize  yields  from  Ar.  Kr,  and  Xe.  The  resulting 
c^Mbiiities  frcmi  low  and  high  coiiq)ressioo  pinches  are  then  conq>ared  with  the  required  source 
yields  in  tabular  form.  The  results  in  the  moderate  to  high  energy  (>5  keV)  photon  range  point 
to  a  fundamental  task  for  load  i^iysics  on  JUPITER:  the  need  to  achieve  ^tenfold  compression 
in  initial  to  final  pinch  radius.  The  last  Secdtm  (VI)  contains  a  summary  of  foe  present  findings, 
liinirariniK  on  tli#»  fvinrliwinnn,  »nA  pmpnsed  directiQna  for  continued  research  on  JUM  i  kk  loads. 


n.  The  2PIMP  Simnlatioo  Code 

ZPIMP  is  a  ^indfically  qrmmetric,  1*0.  radiatkxMnagnelohydrorfynaimc  simulation  code 
for  modeling  z-piiichiiiqilosiooswifo  a  hmqpedcirciiit  model  for  the  generator.  The  simplifying 
feature  of  foe  ZPIMP  code  is  foe  treatment  of  the  main  (dasina  load  as  a  sin^  zone  fying  between 
an  inngr  rsdius  JZi  and  an  outer  radius  Hus  zone  carries  a  mean  ko  (dectroo)  density 

ni(n«),  a  mean  km  (electron)  temperature  a  mean  ionization  level  <  Z  >s  n«/n<,anda 

qwdficknizalkxifoxdtaddh  internal  energy  CiMw  Kgurc  I  shows  the  plasma  nomeoclature  and 
foe  geometric  variables.  The  dynamic  equations  ate  found  by  volume  integration  of  foe  stamlard 
momentum  and  energy  equatioos.  Details  of  die  formulation  can  be  found  in  Rd.6.  Only  a 
summary  of  die  naodel  equatioos  are  presented  here.  The  squrrate  nxmaentum  equations  fm  the 
inner  and  outer  radii  are,  reflectively, 

— (— -j-v®)  =  (pi +p*)2irJlo  +  f  JxS^2irrdr,  (1) 

«  '4  4  c 

and 

—  t;^)  ss  — (pi  +p«  —  Pc)2irJJi  —  ~  T  /  JMS^27rrdr.  (2) 

Ifere  M  is  the  total  shell  mass;  i  the  pinch  length; ««(«,  )  foe  velocity  of  the  outor  (inner)  radios; 

Pi  the  ion  pressure  =  nikaTi  (ka  is  Boltzmann’s  constant);  pe  the  electron  shell  pressure 
=  nJeaTt\ Pc  a  core  pressure  for  die  regimi  between  die  axis  and  iZj;  J*  the  axial  current  density; 
and  .9^  the  •^muthni  magnetic  fidd.  A,  is  the  mic^int  of  the  plasma  shdl.  Because  the  finite 
differencing  is  essentially  Lagrangian.  both  the  plasma  shell  (subscript  a)  and  the  core  region 
(subscript  e)  have  a  tensor  artificial  viscosity  [Ref.7]  of  the  form 

(3)  J£a 

„  tnist  I  Special 


□  □ 


where  is  the  non-dimensional  viscosity  parameter,  p  is  the  mass  density  of  the  zone,  Av 
is  the  velocity  difference  between  the  interfaces  of  the  zone,  and  is  the  zone’s  sound  speed 
=  v/VsCp//*)-  If  the  zone  is  e^qpanding,  9«j«is  set  to  zero.  The  parameter  is  adjusted  in  the 
next  section  to  obtain  ovuall  agreement  with  SATURN  aluminum  data. 

The  equations  fw  die  ion  and  electron  internal  energies  in  the  plasma  shell  are  given  by 

^ =  -J»i2w(i2o««  -  R(Oi)  -  -  Vi)  +  Qieir(I^  ~  i^),  (4) 

and 

^(<Z>  ^  Y««)  =  -p«2ir(J2«i>*- riJ^2irrdr^  (6) 

Thetnmspoctooeffidentsfordiekm-dectroDdiennalequilibntionnteQj.  and  die  resistivily  17 
are  taken  fiom  Braginsldi  [Ref.8].  The  next  to  last  tenn  in  eqn.(S)  is  die  resistive  headng  rate  in 
die  sMI,  and  Krmd  is  the  p]asnia*s  eadiadon  emisskn  rate.  An  ideal  gas  law  is  assumed  and  die 
^pedfie  kmization/excitatioa  energy  eiMu  is  part  of  die  dectron  internal  energy: 

3  » 

ci  ^^haTi/mi;  (6a) 

The  core  r^cm  inade  Rt  is  agMgnftd  1%  of  die  diell  mass,  has  the  same  initial  tenqierature 
as  the  shell,  and  is  treated  as  a  radiadonless  gas.  Its  energy  equation 

d  3 

^(2*’-^Pe)  ~  —pc2irRiVi  ~  (7) 

reflects  bodicompressional  and  viscous  heating.  The  initial  2*IMP  code  described  in  Ref.6  did  not 
include  a  core  zone.  Tmaging  diagnostics  on  a  Russian  alurmnum  wire  experiment  [Ref.9]  indicate 
the  formation  of  a  hot  cote  just  prior  to  plasma  assembly  on  axis.  Detailed  multi-zone  radiation- 
magiMfnhydmrfynafnic  simiilatiftna  of  aigon  inylosions  on  DOUBLE  EAGLE  [Ref.lO]  USedalow 
density  inti#»r  core  piaiema  in  matrJiiiig  die  spectTsl  data  for  the  mean  d^isity  and  tenqierature  at 
aragnarinn.  TtM»  importance  of  tills  ftynamic  feature  was  evidenced  earlv  in  OUT  present  research  bv 
the  failure  to  match  a>IMP  calodations  without  a  cote  to  SATURN  aluminum  data  presented  in 
the  next  section.  Hmce  the  presoit  ZPIMP  version  includes  a  simple,  but  adequate,  treatment  in 
eqn.(7)  of  the  hot  core  formation  in  order  to  better  reflect  the  observed  stagnation  physics. 

The  current  density  and  magnetic  field  ate  related  through  Ampere’s  law 
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The  magnetic  field  obeys  an  induction  equation  fonned  by  combining  a  generalized  Ohm’s  law 
with  Faraday’s  law: 

^  +  W 

In  ZPIMP  these  field  quantittes  can  be  calculated  mi  a  fine  grid  and  fiien  Maxwell  stresses  and  the 
resistive  ue  int^rated  over  different  volumes  as  shown  in  eqns.(l),  (2),  and  (5). 

The  ionization  ^namics  and  radiatian  trauqxxt  is  time  ^lit  from  the  magnetohydrodynamics, 
Le^  die  ^wdficiooizationAexcitationeoeigycjMu*  the  mean  charge  state  <  Z  >.andtheradiati<Hi 
emissioa  rate  Ar«i4  in  die  plasma  shell  zone  are  held  constant  during  a  hydro  timestqp.  At  die  end 
of  a  hydro  dmestqi  and  befoie  die  neact  one,  die  above  three  quantities  are  iqidaied,  subject  to  the 
oooservatioa  of  total  dectroo  intanal  eoeqy  (fib). 

TUpifirfi  impIftdniM  with  fhiirdiffiemnf  eleinmta  wiB  hecooridered:  aluminum,  argon,  krypton, 
and  xentm.  For  die  first  duee  maumml*  die  atomic  level  structure  is  defined  duou^ 

Detailed  Configuratian  Accounting  (DCA).  For  aluminum,  die  bare  nucleus,  all  ground  state  levels 
and  121  eatdted.kvels  spread  tfarous^  the  upper  eleven  ionization  stages  ate  foUowed.  Free- 
bound  transitions  to  and  firam  all  of  diese  levds  are  possible  and  2fi9  bound-bound  transidons  are 
oonddered.  The  imndierB  fior  argon  are:  stripped  4>e^bteen  ground  states  + 108  exdted  levels  in 
the  tqiper  eleven  stages,  and  389  bound-bound  transitioiis;  and  for  krypton:  strqiped  +  diit^-six 
ground  stales -I- filexdled  levels  in  die  upper  deven  stages,  and  184  bound-bound  transidofis.  For 
eadi  s^Miale  material  die  level  populations  are  determined  a  set  of  atomic  rate  equations 

i 

^diere  /u  is  die  finctional  population  of  level  k  in  die  t’th  ionic  species,  and  tftjki  is  the  net 
rate  the  transition  in  die  t’th  species  ffom  the  initial  level  j  to  the  final  level  k. 

The  net  rale  includes  collisimial  ionization,  excitation,  and  de-excitation,  3-body  and  radiative 
fifwmihinatinn,  phoro-ioniTatinn,  pumping,  radiative  decay,  and  inner-sbell  absorption.  In  general 
thft  timn  dfrivati^  mglectod  that  collisional-radiative  equilibrium  is  achieved.  Details  on  die 
solution  tnrhniqiies  and  rate  coefficients  are  summarized  in  Duston,  et  al.  [Ref.l  1]. 

Fbr  xfnmi.  the  dntaili^  configuration  accounting  is  a  potential  future  development  Instead,  an 
Average  AUnn  (AA)  ap{»oadi  is  used  wherein  electrons  take  on  non-integral  occupation  numbers 
within  a  scTpffned  hydrogenic  model  fm  the  atomic  structure.  Our  AA  model  includes  31  levels: 
all  quantum  sublevels  from  n  =  1  to  7,  and  lunqied  sublevels  for  n  =  8, 9,  10.  Thirty-one 
recombiiuuion  edges  and  199  dipole  allowed  tranations  are  possible  with  this  AA  level  structure. 
The  formulas  for  the  electronic  energy  level  structure  are  taken  from  Perrot  [Refs.l2  and  13], 
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with  the  screening  constants  from  More  [Ref.l4]  and  Marchand,  et  al.  [Ref.  15].  The  atomic  rates 
which  couple  different  levels  as  in  eqn.(10)  and  the  radiative  transition  probabilities  are  based  on 
the  hydtogenic  approximation  of  Post,  et  al.  [Ref.l6].  The  present  AA  model  has  been  under 
development  for  sevnal  years  and  represents  a  significant  improvement  over  the  standard  model 
with  only  tm  lumped  tevels. 

Radiation  emission  from  die  plasma  and  opacity  are  dependent  upon  the  local  atomic-level 
populations,  for  not  only  are  photons  created  by  radiative  recombination  and  decay  transitions, 
such  photcms  also  lead  to  population  ledistribudmL  Thus  die  ionization  dynamics  and  radiation 
transput  is  a  strongfy  coupled  problem  and  must  be  solved  togedm:  The  radiation  transport  of  the 
bound-bound  and  bound-ftee  transhioos  are  carried  out  using  the  i»obability-of-escq)e  formalism 
described  in  Aimizese  [Refs.!?  and  18].  Muiti-fiteqneiK^  transport  is  performed  for  die  ficee-fiee 
radiation.  Because  the  tranqiott  {dace  across  only  oae  zone,  the  AA  method  f<v  ionization 
rtynamics  can  use  a  rimilar  radiaticHi  routine  as  the  DCA  method. 

The  generator  driving  the  plasma  load  is  modded  in  ZPIMP  with  a  lumped  circuit  For  the 
ciicuits  in  this  pi^ier  die  wtarhin^.  is  described  an  qpen  circuit  voltage  pulse  Voe*  an  impedance 
and  a  firemt  end  inductance  L«.  From  the  circuit  viewpoint  the  load  omsists  of  a  variable 
inductance  pinch  of  radius  iZe  and  lengdi^  surrounded  by  r^um  current  posts  at  R^: 

+  =  (11) 

In  cgs  units  with  fio  =  4a-/c^,  where  c  is  the  speed  of  light  the  vacuum  magnetic  field  B  is  related 
to  die  load  current  through  B  =  H/rc.  The  initial  load  inductance  is  given  by 

The  generator  is  coupled  implicidy  to  the  plasma  (fynamics  [eqn.(  1 )]  and  the  field  diffusion  equation 
[eqn.(9)].  The  magnetic  and  total  plasma  energy  in  the  load  region  agrees  with  the  time  int^rated 
powo'  delivered  to  the  load  firom  eqiL(l  1)  to  bettn  than  1%. 


6 


Fis.1  Schematic  diagram  for  the  z-pinch  modelmg  in  the  ZPIMP  numerical 
simulation  code.  The  inno'  radios  iZr,  with  velocity  Vi,  and  outer  radius 
Rc,  with  velocity  Vo,  surround  the  iixq>loding  plasma  shell  with  total  mass 
M,  mean  i<Hi(electr(Hi)  drasity  nj(n«),  mean  ion  (dectron)  temperature 
Ti(To),  specific  ionization/excitation  energy  cjonx.  radiative  emission  Arad, 
ariiniithal  rnagnetic  field  B^,  zai  axial  current  density  J^.  The  latter  two 
quantities  can  have  a  gradient  inside  of  Ho.  The  r^on  between  Ho  and  the 
return  cunoit  radius^  is  a  vacuum.  The  radiationless  core  region  between 
Ri  and  the  axis  has  a  pressure  Pe«r«-  Both  the  main  plasma  shell  and  the 
cote  ate  subject  to  an  artificial  viscosity  goi,. 
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m.  Comparison  of  ZPIMP  Calculations  with  SATURN  Af  Data 


A  laige  set  of  aluminum  wire-array  experiments  were  performed  on  SATURN  over  the  past 
few  years.  The  arrays  were  2  cm  long  and  consisted  of  24  wires  of  various  thickness.  The  initial 
load  configurations  ranged  over  an  order  of  magnitude  in  mass  and  nearly  that  in  array  radius. 
The  otperimental  data  including  the  Ai  K-shell  yidd  are  rqxMted  on  and  compared  with  radiation 
scaling  law  jaedictimis  by  Whitney,  et  aL  (Ref.l].  In  this  section  an  analogous  comparison 
is  presented  with  the  ZPIMP  code.  In  terms  of  aun  equivalent  Thevenin  circuit,  the  SATURN 
generalKV  has  a  firemt  end  inductance  of  £•  =  9.75  nH  and  impedance  of  =  0.167(2  (Ref.l9]. 
The  dicuit  and  open  circuit  are  diqdayed  in  Hg.2.  The  load  mass  JIf,  initial  radius  iZ»(i  =  0), 
and  length  i  are  dioseo  to  maldi  die  eaqreiimental  data  set  Tb  start  off  die  present  ZPIMP 
calciilations,  the  initial  load  inductance  of  eqrL(12)  is  fixed  at  2J  nH  1^  adjusting  for  each 
ipedfied  s:  0).  The  inner  radius  .Rr  at  i  =  0  is  always  taken  as  */*Rc{t  ~  0).  The  initial 
tenqMrature  for  die  shdl  and  die  core  plasma  is  S  eV. 

Before  cooqiaring  die  yidd  data  widi  die  calculated  values,  amsider  the  simulation  details  for 
die  iiiqilosioa  widi  M  s  0.657  mg  and  JR«(i  s  0)  >  1.21  oil  Several  properties  of  the  pinch  are 
diqdqred  in  Rgs.3a  »  3d.  The  inductive  notch  of  die  load  current  is  near  die  peak  of  the 

KoaheU  ndiadon  pulse  in  Fig3a,  the  miniriuimof  the  plasma  shdl  radii  in  Fig3b,  the  maximum 
diell  denddes  in  Fig  Jc,  and  the  maximum  diell  temperatures  in  Fig3d.  The  in^losion  time  can 
be  defined  as  the  difference  betwmi  two  fiducial  marks:  (i)  the  intercqit  of  a  linear  fit  to  the  rising 
current  with  the  time  axis,  and  (ii)  the  point  of  minimum  during  the  inductive  notch.  For 
the  case  in  Fig  J,  the  resulting  irrqilosion  time  is  ~S2  nsec.  The  inner  radius  never  touches  the 
central  axis  in  Hg.3b  -  the  whole  plasma  shell  bounces  off  the  hot  irmer  core  during  the  implosion 
{diase.  Note  the  delay  in  terrqierature  equilibration  between  the  ions  and  electrons  in  Fig  Jd.  The 
thermalization  of  the  kinetic  oietgy  in  the  irrqilosion  leads  first  to  a  rapid  and  large  heating  of  the 
ions,  characteristic  of  a  strong  shock  phenomenon.  The  centroid  of  the  K-shell  radiation  pulse 
is  shifted  toward  the  expanding  phase  of  the  pinch,  where  the  electrons  have  come  into  close 
equilibration  with  the  ions.  The  mean  ionization  level  in  Fig3b  rises  to  the  He-like  plateau  at  ^^70 
nsec  and  becomes  nearly  stripped  during  the  highest  conq)ression.  For  this  simulation  the  time 
int^rated  radiation  power  from  tiie  aluminum  K-shell  was  S2.S  kJ.  The  experimental  value  for  the 
same  initial  conditions  was  63.3  kJ. 

Figure  4  displays  the  K-shell  yields,  Yk,  from  the  SATURN  data  set  within  circles  at  the 
q>propriate  location  in  the  Ro  —  M  plane.  The  Yk  value  reflects  the  radiation  measured  through  a 
Vs  mil  Ktqiton  filter.  The  broad  coverage  of  the  SATURN  data  over  the  M  -  Ro  plane  presents  a 
rhallftnging  data  set  to  match.  Unfortunately  none  of  the  shots  were  repeated,  so  the  intrinsic  yield 
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fluctuation  from  shot  to  shot  is  unknown.  In  the  ZPIMP  simulations  the  viscosity  parameter 
of  eqn.(3)  was  adjusted  so  that  the  same  value  provided  a  good  match  to  the  tluee  experimental 
shots  at  Af  =  0.657  mg.  For  the  displayed  fit  in  Fig.4,  =  10.  The  rest  of  the  simulations  were 
performed  then  with  the  same  value.  The  triangles  of  Fig.4  contain  the  K-shell  yields  as  calculated 
from  the  ZPIMP  code  using  the  same  initial  conditions  as  the  experimental  circles  to  which  thqr 
are  attached.  The  K-shell  yields  in  the  simulations  are  the  sum  of  ail  photons  above  1.5  keV,  which 
is  just  below  the  He-pt  line  of  aluminum.  One  immediately  notices  the  remaricably  close  agreement 
between  experimoits  and  simulations  throughout  die  Af  —  .Ro  plane,  except  to  the  right  of  the  plot 
In  this  lattia-  r^on  the  calculations  indicate  that  the  kinetic  energy  of  the  plasma  at  in^losion  is 
insufficient  to  thetmalize  the  entire  plasma  load  into  the  K-shell  ionization  stage.  In  the  parlance  of 
Ref.l,tlieseareimplosi<uiswithi7;S  1.  Tlielargeroqierimental  values  in  this  domain  compared  to 
those  from  simulations  arise  from  the  presence  of  large  qiatial  gradients  in  the  plasma  at  inqilosion. 


Zos0.167a  Lo*s9.75nH 


SATURN 


Fig.7  Equivalent  Thevenin  circuit  for  the  SATURN  generator  and  the  open 
circuit  voltage  K»e  waveform.  Fbr  the  calculations  of  the  present  paper  the 
initial  load  inductance  is  fixed  at  2.5  nH. 


9 


:  (MA) 


t  (nsec)  t  (nsec) 


Fig3  Dynamic  features  of  an  aluminum  implosion  on  SATURN  from  the 
ZPIMP  simulation  code.  Initial  conditicms  are:  M  =0.657  mg,  initial  .Ro  = 
1.21  cm,  initial  Ri  =  *URo  cm,  and  pinch  length  /  =  2  cm.  The  viscosity 
parameter  is  =  10.  (a)  load  current  and  K-shell  radiated  power,  (b) 
outer  and  inner  radius,  with  Vio  the  mean  ionization  level  <  Z  >',  (c) 
electron  and  ion  densities;  (d)  electron  and  ion  temperatures.  The  calculated 
K-shell  emission  is  52.5  U  while  the  experimental  value  for  the  same  initial 
load  configuration  was  63.3  kJ. 
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initial  radius  (cm) 


mass  loading  (mg) 


Fig4  Aluminum  K-shcll  yields  for  SATURN  in  the  configuration  plane  of 
load  mfliK  -  initial  radius.  The  dicles  contain  the  expeiiniental  data  in  kJ 
from  2  cm  l<mg  wire  arrays.  The  triangles  contain  the  K-shell  yields  as 
calculated  firom  the  ZPIMP  code  using  the  same  initial  conditions  as  the 
experimental  circles  to  which  they  ate  attached.  The  viscosity  parameter 
10. 
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IV.  Krypton  Radiation  Yields  on  a  JUPITER-Class  Generator 


In  the  present  paper  an  early  design  version  of  the  linear  inductive  voltage  adder  [Ref.20]  is 
usedastheJUPirER-dassgeDerativ.  A  transmission  line  code  was  used  to  convert  the  design  into 
an  equivalent  Tbevenin  lun^wd  element  drcuit  The  resulting  parameters  and  q)en  circuit  voltage 
Voc  for  diis  JUPIIER  generator  is  displayed  in  Fig.5.  To  stutty  the  dqiendenoe  of  the  yield  over 
a  range  of  peak  load  currents  we  scale  the  machine  throu^  the  opca  circuit  voltage,  kequng  the 
odier  circuit  parameters  fixed.  Spedfically,  the  actual  open  circuit  voltage  is  given  by  aVmct  with 
0.4  <  tt  <  2.  The  cooesponding  machine  will  be  denoted  by  LIA-l(a).  Hereinafter;  the  2PIMP 
calodatkMis  for  JUFTTEk  assume  /  s  4  cm.  =  2.63  nH.  »  0)  s  =  0).  and 

Ti{t  =  o)  *  T.(l  *  0)  *  6  eV 

Tn  nwWtft  ftotermingthe  Inad  maan  and  ttiitial  r^ingeftndifinn*  fn  npd"**!  Irtypr^ 

radiatioa  yields  it  is  not  practical  to  perform  hundreds  of  ZFIMP  runs,  gradually  narrowing  down 
the  optimal  r^<m  of  parameter  qrace.  Instead,  a  radiation  scaling  law  is  used  to  pve  a  broad  brush 
incture  die  K-shell  yield  dependence  in  the  Jlf  —  A0  plane.  A  diin  shell  model  fm  the  fnndi 
dynamics  rqrlaoesZPlMP^  die  impkadon  is  stopped  once  die  ciicuit  has  driven  dreshdl  to  ^/r 
ofits  initial  radius.  Details  ofdie  procedure  are  described  in  ReO.  Krypton  K-shellyidd  contours 
fordie  JUPITER  generator  LIA-l(a  s  l)aredi^dsyedinRg.6.  ItsuggestsdiataseriesofZPlMP 
runs  to  obtain  irrquroved  yidd  esdirutes  should  be  performed  with  initial  radii  firom  3  tti  7  cm  and 
load  masses  firom  io  to  40  mg. 

Ofthe  six  amuladons  performed  die  peak  kryiXonK-shell  yield  for  LIA-l(a  1),  only  SOU, 
isfoundat  JIf  ^  20ingandHo  Scol  Consider  some  details  of  the  dynamics  for  diis  simulation. 
The  value  for  the  viscosity  parameter  is  the  same  as  that  used  to  match  the  SATURN  data 
(slO).  Fig.7  shows  die  evolution  of  die  load  current,  K-shell  radiated  power,  iruier  and  outer 
radii,  mean  itmizadon,  imi  and  electron  density,  and  ion  and  dectron  teneperature.  l^asuring  the 
inplosion  time  in  a  similar  way  as  diat  described  for  SATURN,  Fig.7a  indicates  that  Ump  ~  130 
nsec.  One  can  see  (Fig.7a)  that  the  inductive  notch  of  the  load  currrat  has  dropped  the  current 
from  a  maximum  of  ~SS  MA  to  ~33  MA,  die  outer  radius  of  the  plasma  load  has  decreased  only 
by  afactorof 'wS,  fromRo(f  -  0)  =  5  cm  to  1  cm  (Fig.7b).  and  the  peak  dectron 

tmrperature  (Hg.7d)  is  ~3  keV  at  implosion.  The  ion  temperature  spike  in  Fig.7d  represoits  the 
thennalization  of  die  Idnetic  enogy  at  implosion. 

Figure  8  shows  the  2PIMP  time  int^rated  tynthetic  spectra  for  the  above  simulation  run.  The 
forest  of  narrow  emission  lines  below  5  keV  arises  from  L-sbell  bound-bound  transitions.  The 
recombination  edge  to  the  ground  state  of  Li-like  krypton  lies  between  these  lines  and  the  5  keV 
marie.  The  strong  He-a  line,  a  ls2p(^Pi)  to  ls’(*So)  transition,  is  clearly  seen  at  <>.13  keV.  This 
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energy  is  the  lower  end  of  the  K-shell  radiation.  No  other  prominent  K-shell  emission  lines  are 
seen  for  this  case.  The  continuum  beyond  '^'17  keV  is  the  sum  of  free-bound  transitions  to  the  He- 
and  H-like  ground  state.  As  mentioned  above,  the  K-shell  yield,  i.e.,  the  integrated  emission  >  13 
keV,  for  this  simulation  is  50  kJ.  Note  that  the  He-a  line  is  a  major  contributor  to  this  yield. 

In  Fig.9  the  yields  within  the  various  energy  bins  marked  on  Fig.8  are  shown  above  the  SS 
MA  peak  load  current  point  The  yield  in  the  15  -  20  keV  bin  is  included  in  the  K-shell  yield, 
but  <mly  ~  of  the  yield  in  the  10- 15  keV  bin  crmtributes.  lb  di^lay  the  trend  of  yield  with 
load  current,  a  similar  search  as  done  above  was  performed  for  LIA-l(a  =  ^^)andllA-l(a  =  2) 
circuits.  For  each  of  these  madiines,  die  scaling  law  was  first  used  as  in  Hg.6  to  provide  a  tough 
picture  of  the  yield  dqiendenoe  <m  die  initial  load  configurafimis.  Frmn  diese  |dots  a  limited  set  of 
M  —  iZ0(<  as  0)  conditions  were  chosen  for  fiirther  sturty  widi  ZPIMP  simul^ons.  The  rqrtimal 
conditions  were  found  to  be  JIf  =  4  mg.  Ro{i  ss  0)  =  3  cmforLlA-l(a  ss  0.4),  and  JIf  =  40  mg, 
RJ^i  =  0)  =  6  cm  for  LIA-l(a  s  2).  The  calculated  yidds  in  the  same  {rfioton  enet^  ins  for 
diese  two  runs  are  also  presented  in  Rg.9. 

For  a  JUPITER  driver  of  ~60  MA  peak  load  current,  one  is  hoping  to  achieve  much  hi^ier 
yields  in  the  10-15kBVbinthan8howntnHgil.  Ibptoduce  1  Ml  <^10-15  keV  photons  this 
figure  suggests  dial  it  is  necessary  to  devdt^  a  lUPTIER  machine  wfaidi  can  ddivei  at  least  80 
MA  pe^  load  current,  a  formidable  task. 

The  essential  problem  in  {xoducing  hitter  radiation  yidds  is  the  lack  of  adequate  conqnessimi 
at  inqilosion.  The  compression  ratios  of  Roi^  —  0)/Ao(t{mr)  ^  above  JUPITER  calculations 
aresiniilartothe  Vs  value  found  in  the  SATURN  simulations.  The  choice  ofthe  viscosity  parameter 

which  essentially  controls  the  conqrtession  ratio,  was  based  on  matching  3*IMP  with  the 
McperiitiMital  aluminum  K-shell  yidds  for  SATURN.  Suppose  it  is  possible  to  achieve  higher 
compressions  on  JUPITER  than  on  SATURN;  how  would  that  affect  the  yidds?  To  study  this 
question  requires  a  hardo'  inqilosion  in  die  simulations,  ix.,  one  in  which  the  thmnalization  of 
enet^  is  delayed  until  the  pinch  is  tights.  This  means,  as  far  as  the  simulation  modeling 
is  concerned,  that  the  viscosity  needs  to  be  lowered.  The  same  M  and  Ro[t  =  0)  initial  conditions 
listed  in  Fig.9  were  repeated  but  widi  the  viscodty  parameter reduced  to  3  instead  of  the  value 
10  used  above.  The  results  for  the  yield  variation  with  load  current  is  presented  in  Fig.l0.  An 
aAiitiftnai  point  at  ~40  MA  peak  load  current  is  also  included.  The  yields  in  high  energy  bins 
above  10  keV  have  dramatically  increased.  To  reach  1  MJof  10-  15keV  photons  onUA-lfa  =  1) 
now  can  be  achieved  with  55  MA  peak  load  current,  sigiuficandy  less  than  the  '>^80  MA  found  in 
the  low  compression  case  of  Fig.9.  The  K-shell  yields  in  this  high  compression  scenario  are  much 
closw  to  the  iyaling  law  predictions,  as  contained  for  instance  in  Fig.6.  The  K-shell  yield  for  the 
M  =  20  mg,  Ro{i  =  0)  =  5  cm  tun  is  now  1040  kJ. 
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To  discuss  the  stagnation  physics  which  is  leading  to  the  improved  yields,  consider  the  details 
of  the  high  compression  implosion  in  Fig.  1 1.  The  initial  conditions  (M  =  20  mg,  iZe  =  5(<  =  0) 
cm)  are  the  same  as  shown  for  the  low  compression  case  shown  in  Fig.7.  Though  the  same 
peak  load  cuxroit  is  reached,  note  the  deeper  inductive  notch  in  Fig.  1  la  and  the  associated  (toeper 
in^losion  in  Hg.llb.  Because  the  viscosity  is  lower  than  in  Fig.7,  the  bounce  is  harder  and  the 
peak  kinetic  energy  is  larger  ~11  MI  in  die  hi^  con^nession  case  conqiared  to  ~8  MJ  in  the 
low  one.  In  both  cases  the  ion  teoqietature  shoots  up  to  tens  of  keV  at  thennalizatimi,  but  the  peak 
electrcm  ten^ierature  is  ~S  ke V  for  hi^  amqxession  case  (Fig.  lid),  and  only  ~3  ke V  for  the  odier 
(Hg.7d).  Furthermore,  die  electron  densities  are  in  die  rado  ne(bigb)^Gow)  =  Figs.7c 

and  11c.  Ibgetber.diesefBcdKS  lead  to  a  shorter  ion-electnmequilibradoo  time. 


TU  oc 


T« 


<Z>n^* 


(13) 


for  the  high  conqitessioa  case.  This  is  consistent  with  die  more  rqnd  equilibration  and  hi^ier  peak 
T«  evidenced  in  Fig.7d. 

The  inqnct  upon  die  K-shell  yield  due  to  the  larger  fic  and  Tc  for  die  hi^conqiression  case 
can  be  seen  by  noting  die  formula  fm  the  He^  emission  rate: 


sec 


g-is(*«v)/r. 


(M) 


whme  liHe  is  the  He-like  ground  state  density  of  krypton.  This  equation  is  valid  in  the  limit  of 
corcmal  equilibrium  where  evoy  upward  excitation  to  the  excited  state  is  followed  by  a  radiative 
decry.  To  attain  die  opposite  extreme  of  local  thermodynamic  equilibrium  in  the  K>shell  of  krypton 
requiresnc  ~  10**  cm~*,whidt  far  exceeds  the  peak  ne  in  Fig.llc.  This  large  value  follows  from 
the  correspoMiding  required  electron  density  for  aluminum  (10**)  [Ref.21]  and  the  ZJ^^  scaling  to 
other  elements  [Ref. 22] .  Using  the  above  values  for  Te  and  leads  to  a  factor  of  ~20  enhancen^nt 
in  the  He-a  bound-bound  photon  production  rate  fm  the  high  conqiression  case  compared  to  that 
for  low  conqnessimL  This  factm  would  be  evra  largm'  if  the  slight  increase  in  mean  ionization  of 
Fig.  1 1  b  over  that  in  Fig.7b  were  taken  to  indicate  a  slighdy  enhanced  n^c  ftn*  the  Mgh  coni(nession 
case.  Given  that  the  K-shell  pulse  widths  in  Fig.7a  and  11a  are  nearly  the  same,  and  the  dominance 
of  the  He-a  contributitm  to  the  total  K-shell  yield  for  this  implosion,  the  factor  of  ~20  explains  the 
ratio  in  the  K-sbell  yield  of  the  two  cases:  yjr(high)/yjir(low)  =  (1040  kJ)/(50  kJ). 

Thus  the  large  diffmence  in  the  yield  is  due  to  moderate  changes  in  the  plasma  conditions 
during  the  thennalization  phase,  i.e.,  the  faster  temperature  equilibration  coupled  with  the  enhanced 
electron  density  of  a  higher  conqiression.  As  pointed  out  in  the  first  prqier  in  this  series  on  scaling 
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law  estimates  of  JUPITER  yields  [Ref.3].  the  optimal  yielding  implosions  for  krypton  on  JUPITER 
fne  weak  implosions  in  the  sense  that  the  kinetic  energy  is  barely  sufficient  to  thennali7/>  the  whole 
plasma  into  the  K-shell  ionization  stage.  Using  ^  as  the  ratio  of  the  kinetic  energy  per  particle 
over  die  energy  needed  to  ionize  the  atom  to  a  50%He-  SO%H-like  configuratitm  [Ref.4],  these 
inoplosions  have  17  ;$  1.  Such  implosions  were  only  crwiely  treated  by  the  J-scaling  law  described 
in  Ref.3,  and  the  present  results  of  Figs.9  and  10  attest  to  the  sensitivity  of  the  yields  arising  from 
die  complexity  of  the  stagnadtm  physics. 


Zo-O.24a  LoiJ503nH 
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VOLTAGE  ADDER 
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0  100  200  300 
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Fig3  Equivalent  Thevenin  circuit  fm  a  JUPITER  generator.  The  circuit  is 
based  on  an  early  versimi  of  the  linear  inductive  voltage  adder,  herrinafta 
called  UA-l.  For  the  calculations  of  the  present  pspo*  the  initial  load 
inductance  is  fixed  at  2.63  nH  and  the  pinch  length  is  4  cm.  For  differmt 
machine  versions  the  anplitude  of  the  open  circuit  voltage  is  multiplied  by 
a  parameter  a  between  0.4  and  2.  The  corresponding  version  is  termed 
LIA-l(a). 
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E-Shell  Yield 

T  iJSUl 

I  1.30B403 


2.1SB+01 


>41 


4.85B««0 
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Flf^S  Yield  oootbuts  in  tbe  load  ni&ss  initial  radius  plane  for  krypton  K-sbeU 
radiation.  Tbe  yields  are  based  on  a  ctnnbinadoD  of  die  thin  Adlttynainics 
forthe  load,  the  LIA-l(a  ==  1)  design  ofFig.5,  and  the  K-sheUJ-scaling  law 
of  Ref3.  The  inqiloaon  is  stq^ied  at  Vr  of  the  initial  racUus.  This  grqih 
suggests  die  r^on  of  initial  conditions  where  ZITMP  calculations  can  focus 
<ni  inqiroved  yield  estimates. 


LIA-Ua  s  l)firoiii the  ZPIMPsiniulation  code.  Initial  ccHiditi<His  are:  total 


mass  M  s  20  mg,  initial  JZo  =  5  cm,  initial  Ri  —  */4Ro  cm,  and  pinch 
lengdi  i  —  4  cm.  The  viscosity  parameter  10)  is  chosen  to  give  a 
similar  con^xession  (Vs)  as  used  to  match  the  SATURN  data  in  Fig.4.  (a) 
lottl  current  and  K-shell  radiated  power;  (b)  outer  and  inner  radius,  with 
Vs  of  the  mean  ionization  level  <  Z  >;  (c)  electron  and  ion  densities;  (d) 
election  and  ion  tenqieratures.  The  calculated  K-shell  emission  is  SO  U. 
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Intensity  (ergt/cnrke?) 


peak  load  current  (MA) 


Elf.9  Kiyplon  yidds  vs  peak  load  current  in  various  photcm  energy  bins  for 
die  dicahs  LIA-l(o  »  0.4,1.0,2.0).  The  load  mass  and  initial  radius  f<v 
the  differentiii9losi(»s  are  shown  the  top.  The  viscosity  parameter 
10)  for  this  case  leads  to  low  compressions  where  i2o(f«mp)/<Ro(t  =  0)  is 


yield  (kJ) 


2mg  13  20  40 
3cm  4  5  6 


peak  load  current  (MA) 


F%.10  Kiypton  yidds  vs  peak  load  cunrat  in  various  photon  energy  bins 
for  die  dicoits  UA-Ua  —  0.4, 1.0, 2.0).  The  load  mass  and  initial  radius 
for  the  different  in^ylosions  are  shown  at  the  top.  The  viscosity  parameter 
3)  for  diis  case  leads  to  high  conqnessionswhaeJ2o(<«mp/-Ra(<  =  0) 
is  ~  */io. 
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0  100  200 


Flg.ll  Dynamic  features  of  a  krypton  inqplosion  on  the  JUPITER  ciicuit 
LIA>l(a  s=  1)  from  the  ZPIMP  simulation  code.  Initial  load  conditions  are 
the  same  as  in  Fig.7.  The  visc(»ity  parameter  3)  is  diosen  to  give 
a  higher  ctmqnession  (~  Vio)  than  in  Hg.7.  (a)  load  current )  and  K-shell 
radiated  power;  (b)  outer  and  inner  radius,  with  Vs  niean  ionization  level 
<  Z  >;  (c)  electron  and  ion  densities;  (d)  electron  and  ion  temperatures. 
The  calculated  K-shell  emission  is  1040  kJ. 


21 


V.  JUPITER  CqMbilities  vs.  RadiatioD  Requirements 


The  final  task  before  us  is  to  use  the  ZPIMP  simulation  code  to  predict  JUPITER  yields  in 
various  energy  bins  and  compare  the  results  with  stated  radiation  testing  requirements.  Table  I  lists 
radiation  fluence  requirements  for  simulator  testing  and  the  corre^ndmg  source  yields  resulting 
from  fluenoe  unif<»mity  on  target  The  inqiact  on  source  yield  due  to  various  standoff  distances 
inq)osed  debris  shielding  is  presented  in  Ihble  IL  Thble  m  lists  the  optimal  yields  in  the  energy 
bins  from  argtm,  krypton,  and  xenon  calculated  with  the  ZPIMP  code  for  the  JUPITER-class 
machine  LIA-l(a  1)  of  Fig.S.  A  detailed  desct^tm  of  the  tabular  columns  foUows. 

FrmnT.  Kermetty's  presentation  at  JDOST-3  ptefl23]  die  fluence  requirements  of  Thble  I  are 
sqiantfed  into  five  (dioton  energy  bins,  denoted  by  die  letters  A  -  E.  The  bin  nomenclature  and 
energy  limits  are  listed  in  die  first  two  columns.  Cases  D  and  E  have  been  amended  to  include  an 
iqiper  bound  on  the  photon  energy.  The  minimal  requirements  from  T.  Kennetty  of  fluence  on  die 
target  surface,  target  area,  and  uniformity  are  given  in  the  third  dirough  fifth  column,  reflectively. 
Assume  that  die  target  presents  a  disc  normal  to  die  line  of  sight  from  the  source.  The  radius  of 
such  a  disc,  r,  is  readily  oonqrated  in  column  six  from  die  stated  area.  Let  D  be  the  distance  from 
the  source  to  center  of  die  tuget, «  die  uniformity  requirement,  and  <  the  filter  transmission  factor. 


The  geometry  is  dqikted  in  Fig.12.  If  y  is  die  srxirce  radiation  yield,  the  fluenoe  at  the  center  of 


dietai]^is 


^crater 


Y 


(15) 


vdiile  the  fluenoe  at  the  target  edge  is 


4ir(I?’  +  r*) 


COB  9. 


(16) 


The  angle  0  mises  from  the  fact  that  the  radiatimi  from  the  source  striking  the  edge  of  the  target  is 
at  an  angle  to  its  surface  normal.  Since  the  tat]^  was  assumed  flat,  the  angle  is  the  same  as  that 
subtoided  by  the  target  at  the  source.  Because  the  target  shfies  may  not  be  flat,  results  in  Table 
I  are  given  with  and  witlxnit  the  angle  ocarectioiL  To  satisfy  the  uniformity  requirement  one  must 
have  Fu§*  >  anti  th®  minimum  distance  which  can  satisfy  this  requirement,  n^ecting 

the  angle  term  in  eqn.(16),  is 

■Dmin  =  ^  (17) 

On  the  other  hand,  if  the  angle  term  is  included  and  one  notes  that  coa  6  =  Djy/D^  +r^,then 


Ttw  resulting  D^tn  are  listed  in  columns  seven  and  eight  and  show  that  the  angle  factor  can  make  a 
significant  difference.  Given  the  transmission  factor  of  column  nine,  and  the  distances  of  columns 
seven  and  eight,  the  minimum  required  source  yield  in  each  eneigy  bin  is  presented  in  the  last 
column  of  Table  I.  For  comparison,  we  note  that  the  largest  K-shell  yields  to  date,  76  kJ  above  1 
keV,  has  been  adiieved  on  SATURN  with  A1  wire  array  inq[>losions.  An  inqxovement  of  a  factor 
of  *^22  is  needed  by  JUPTIER  just  to  meet  the  minimum  yield  requirement  fOTcase  B. 

The  minimimi  distance  for  uniformity  of  IU>le  I  does  not,  however,  account  for  the  physical 
iq)paratus  needed  to  protect  the  target  fimn  particulate  debris.  This  apparatus,  wifii  its  shields  and 
fast  valves,  inqioses  a  fd^acal  squration  between  die  source  and  tihe  test  object  vdudi  is  termed 
the  drtnas  mitigation  distance.  The  mitigation  q)patattts  for  a  JUPITER  class  generator  demands 
a  rignificant  R  &  D  effort  Anticipating  future  techndqgy  inqxovements,  ThUe  n  includes  an 
estimated  range  for  dte  debris  mitigation  distance  in  oohunn  five.  The  five  test  cases  A  •£,  tiidr 
energy  bins,  required  fluences,  and  transmissi<»  factors  are  also  rqieated  in  tile  first  four  columns  of 

Tkble  IT.  Baaed  ow  the  flngnee  and  frMinnisritwi  faetor,  and  tiwt  min/maT  rangft  in  AJirig  miripifinn 

distance,  tile  required  source  yields  can  readily  be  calculated  finomeqn.(15).  These  are  presented 
as  rsmges  in  tile  oorreqxndipg  bins  in  the  last  column  of  Ihble  IL  Clearly  tiiere  is  a  lot  to  be  gained 
in  terms  of  the  required  yiehls  if  tile  R  &  D  on  ddifis  mitigation  can  reduce  tile  standoff  distance. 
Note  tl^  the  minimum  debris  mitigation  distance  for  smne  of  tiie  cases  A  •  E  is  smaller  than  tiie 
cotre^Kinding  minimiim  distance  for  unifarmity  based  on  eqn.(18)  and  listed  in  IhUe  L  On  tiie 
odier  hand,  if  tiie  ang^  between  tiie  radiation  and  tiie  target  normal  is  unimportant  tiie  uniformity 
condition  is  autmnatically  satisfied  by  conqiaring  the  distances  of  column  sevm  of  Thble  I  and 
those  in  colunm  five  of  Ihble  IL  A  further  discussion  of  this  issue  is  given  Hettomann  in  the 
Final  JDOST  Repmt  [Ref.24]. 

Given  tiie  source  requirements  of  Ihble  I  or  Tkble  n,  v^iat  can  one  expect  in  tiie  various 
energy  bins  ftmn  a  JUPTIER-class  generator?  An  answCT  to  this  question  can  be  provided  by  using 
ZPIMP  to  calculate  the  radiation  yields  from  specific  load  configurations.  The  specific  JUPITER 
design  used  in  tiie  calculations  is  that  of  Fig.5,  with  a  =  1.  As  above,  a  goieral  view  of  the  yield 
dqpendence  on  initial  radius  iZo(f  =  OjandloadmassJlf  was  created  by  using  the  K-shell  J-scaling 
law  [Ref3]  on  tiiis  generator  fcx’  argcm,  krypton,  and  xenon.  Once  the  promising  region  in  the 
M  —  Ro  was  determined  (the  pinch  loigth  was  fixed  at  4  cm),  ZPIMP  runs  were  performed  for  a 
limifjtH  set  of  ctuiditions  within  this  r^on  and  the  optimal  configuration  was  sought.  For  argon 
and  kryptmi  the  icmization  dynamics  was  treated  with  a  DCA  model  as  described  in  Section  11; 
for  xentm,  the  AA  model.  Because  the  compression  of  the  pinch  based  on  the  viscosity  parameter 
0^  was  found  to  be  very  important  in  the  krypton  yields  of  Section  IV,  both  low  compression 
09^  s  10)  and  high  compression  09«,'«  =  3)  runs  were  made.  The  optimal  material,  load 
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configurations,  and  calculated  yields  are  presented  for  each  radiation  bin  in  Table  III.  Note  the  large 
yield  improvements  with  the  high  (column  eight)  over  the  low  (column  seven)  compression  cases, 
especially  for  bins  C  through  E.  As  explained  in  the  previous  section,  this  difference  reflects  the 
higher  electron  densities,  faster  thermal  equilibration,  and  enhanced  K-shell  emission  of  the  high 
conqnession  runs  coiiq)ated  with  the  low. 

The  final  two  columns  of  Table  01  repeat  the  range  of  required  yields:  column  nine  from  die 
last  column  of  Table  I;  column  ten  from  tbe  last  column  of  Table  n.  The  values  in  column  nine  are 
based  on  eqn.(  1 8),  ie.,  (be  cos  0  term  of  eqtL(  16)  is  included.  Tbe  maximum  in  the  projected  range 
of  required  yields  (columns  nine  and  ten)  is  clearly  that  arising  from  tbe  largest  ddiris  mitigatiem 
distance  -  die  upper  value  of  the  last  column  in  Table  m.  The  minimum  in  tbe  same  jnojected 
rrmge  is  firom  column  nine  for  bins  A  -  C,  and  tbe  lower  value  of  the  last  column  for  bins  D  and 
£.  The  last  four  columns  of  Table  m  offer  a  direct  conqiatison  between  [uedicted  capabilities  and 
anticipated  requirements  for  the  source  yield  on  JUPITER.  For  bin  A.  even  the  simulation  with 
low  compression  exceeds  the  largest  yield  requirmnent  The  radiadon  in  this  lowest  energy  photon 
bin  is  entirely  from  argon  L>shell  continuum  and  line  emission.  The  plasma  Te  need  only  be  a 
coi^le  hundred  eV  to  radiate effidendy  in  die  0- 1  keV  r^ime,  Le.,  the  photon  ouqiut  is  not  ded 
to  a  large  inqilosion  vdodiy.  Ifence  a  large  load  mass  can  be  used  to  boost  the  yield.  FbtbinB. 
only  the  required  yidd  for  the  largest  debris  midgadon  distance  is  Ix^cmd  the  predicted  krypton 
opabilides  for  low  or  higb  oompeesdon  The  next  bin,  C,  presents  a  transidonal  dtuadon:  The 
low  conqnession  Idypton  yidds  ate  bdow  all  of  the  teqmred  yidds,  while  the  opposite  holds  for 
the  high  compression  case.  For  the  high  energy  bins  D  and  E,  both  the  low  and  high  conpression 
simuladons  for  krypton  and  xenon  give  yidds  below  the  range  of  required  values,  though  the  high 
compression  krypton  yidd  for  bin  D,  280  kJ,  is  cl(»e  to  the  required  yield  if  the  minimum  debris 
midgadon  distance  is  realizable. 
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D 


Fig.  12  Geometiy  used  to  transpose  a  radiation  fluence  requirement  at  a  taiget 
of  area  A  into  an  yield  requirement  Y  from  a  point  source  a  distance  D  away. 


Table  I.  Minimum  Yields  for  Uniformity 


1 

ptwton 

bin 

(keV) 

refilled 

floenoe 

(cal/bm*) 

required 

test 

area 

(cm*) 

required 

fluenoe 

unifonniQr 

test 

radius 

(cm) 

minimum 

distance 

leqn.(17)l 

(cm) 

minimum 

distance 

[eqn.(18)] 

(cm) 

filter 

transmission 

factor 

required 
source  yield 
(unifonnity) 
(kl) 

A 

0-1 

1 

80% 

89^ 

178.4 

2093-3261 

B 

1-5 

1 

80% 

79.8 

159.6 

1675-2612 

C 

5-15 

1 

80% 

56.4 

112.8 

0.8 

837-1304 

15-30 

5 

90% 

4.4 

131 

163 

0.8 

57-87 

BSE] 

5 

25 

90% 

18 

8.4 

10.4 

03 
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Table  n.  Mintmum  '^elds  for  Debris  Mitigation 


photon 

bin.. 

OmV) 

- j 

required 

flueiioe 

(cal^*) 

fUter 

traMwii—inw 

fKtor 

min— max 
debris 
mitigation 
distance 
(cm) 

^pied 

somceyidkl 

(ddxis) 

(kJ) 

A 

0-1 

1 

MSM 

200  -  300 

2630-5918 

B 

1-5 

1 

■9 

180  -  270 

2130-4793 

C 

5-15 

1 

125-190 

1027-2374 

□ 

15-30 

5 

0.8 

30-50 

296  -  822 

30-60 

5 

0.8 

20-40 

132  -  526 

Table  m.  ZPIMP  Voltage  Adder  Simulations  [L1A-I(a  =  1)] 


I^iotmi 

bin 

(keV) 

material 

initial 

radins 

(cm) 

pindi 

length 

(cm) 

low 

com{»essi<m 
source  yield 
(kJ) 

high 

compression 
source  yield 
(kJ) 

required 
source  yield 
(unifoimity) 
(kl) 

required 
source  yield 
(debris) 

(kl) 

A 

0-1 

Ar 

60 

3 

a 

7820 

11000 

3261 

2630  -  5918 

B 

1-5 

At 

20 

3 

3060 

4390 

2612 

2130  -  4793 

C 

5-15 

Kr 

5 

609 

2570 

1304 

1027-2374 

D 

15-30 

Kr 

10 

5 

17 

280 

87 

296-822 

E 

30-60 

Xe 

5 

5 

n 

'  - 

50 

36 

132-526 
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VI.  Summaiy  and  Conclusions 


In  a  previous  paper  on  the  JUPITER  PRS  scoping  study  we  noted  three  fundamental  load 
issues:  stagnation  physics,  high  Z  ionization  dynamics,  and  stability.  These  rendered  the  predicted 
K-shell  3tields  from  the  J-scaling  Isw  uncertain.  The  present  report  was  focused  on  providing  an 
in^Koved  estimate  for  yields  in  specific  energy  bios  using  a  1  >D  radiation-magnetohydrodynamic 

code,  2P1MP,  to  simulate  implosions  on  one  JUPTIER  circuit  The  ZPIMP  code  addressed  the 
first  two  of  the  above  issues:  stagnatka  {diysics  by  following  the  load  dynamics  frmn  some  early 
ionized  {diase  tiuougfa  the  in^doaon  and  bounce,  and  ionization/exdtation  dynamics  by  using  a 
detailed  atomic  physics  padcage  witii  sdf-conasteat  radiation  transport  This  model  was  ooiqried 
to  a  sinqdelunqiedcircnit  for  tile  JUPITER  driver.  The  advantage  of  ZPIMP  is  its  sinqde  treatment 
tile  {dasma  load  as  a  single  zone  or  shelL  There  are  sqiatatemoinenmm  equations  fertile  inner 
and  outer  shell  radii,  a  single  ion  internal  energy  equatum,  and  a  single  electron  internal  energy 
equation  (Fig.l).  Since  the  compiitati<»  time  for  radiation  transport  scales  faster  than  the  square 
of  tile  number  of  zones,  this  single  zone  qiptoach  offers  expediency  in  calculatitxL  In  our  study 
of  SATURN  inqdosions  it  was  found  necessary  to  include  a  non-radiating,  low  density,  inner  zone 
which  developed  into  a  hot  core  at  iinplosion. 

T&  vtabili^  of  ZPIMP  was  investigated  1^  comparing  the  aluminum  K-shell  yields  from 
Z'pindi  eaqieriments  on  SATURN  witii  STIMP  simulations.  The  results,  in  tiie  initial  load 
configuration  (Af  —  iZo)  plane  of  Fig.4,  showed  excellent  agreement  witii  the  data,  excqit  in  die 
r^onrrfweakinploskMis  where  the  kinetic  energy  is  insufficient  to  ionize  the  entire  plasma  into  the 
K-shell.  The  viscosi^  parameter  of  the  simulations  (  =  10)  was  determined  by  matching  three 
of  the  data  points  and  then  using  the  same  value  for  the  remainder  of  the  simulaticms.  A  detailed 
stutfy  of  one  of  tiie  implosions  in  Fig  J  indicates  a  tyjncal  conpression  ratio  Ro{Ump)IRo{i  —  0) 
of~  Y*. 

In  cudo’  to  calculate  radiation  yields  from  a  JUPlTER-class  generator,  an  early  version  of 
the  linear  inductive  voltage  adder  (LIA-1)  was  converted  to  an  equivalratlhevenin  circuit  (FigJ). 
The  circuit  was  used  to  rpresent  more  and  less  energetic  machines  scaling  the  <pen  circuit 
voltage  (s  aVoc)  while  keeping  the  otiier  circuit  parameters  fixed.  This  enabled  us  to  look  at 
tiie  yield  dependrace  as  a  function  of  the  peak  load  current  For  eadi  machine  version,  turned 
LIA-l(a),  a  search  was  performed  in  the  JIf  —  iZo  plane  to  find  optimal  yielding  conditions.  A 
typical  implosion  dynamics  was  presented  in  Hg.7.  The  time  integrated  spectra  (Fig.8)  was  divided 
into  various  energy  bins  and  the  yield  in  each  bin  was  plotted  in  Fig.9,  along  with  the  optimal  yields 
for  the  other  machine  versions.  This  figure  showed  disppointing  yields,  and  to  improve  them  the 

parametu,  initially  set  to  match  the  SATURN  data  (=10),  was  lowered  to  3  to  obtain  higher 
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compressions.  The  high  compression  yields  in  Fig.lO  were  significantly  improved.  Comparison 
of  the  dynamics  for  the  high  compression  case  (Fig.ll)  with  the  same  load  conditions  Init  low 
corrqrtession  (Fig.7),  demonstrated  that  the  yield  improvement  was  due  to  a  speedier  equilibration 
of  the  electron  temperature  in  conjunction  with  a  higher  electron  density  at  implosion.  The  moderate 
changes  led  to  a  fadxxof  ~20  enhancement  of  the  total  K-shell  yield,  most  of  it  coming  fipom  die 
He-a  line  emission.  This  result  pointed  out  die  significance  of  stagnation  physics  modeling  for 
those  in^losituis  widi  ^  1 

Finally,  the  2P1MP  code  was  used  to  estimate  optimal  3rields  in  the  energy  bins  listed 
in  die  requirements  for  JUFTIER.  This  seaidi  was  performed  for  LIA-l(a  s  1)  of  Fig.5 
over  the  pmch  tnutwriaU  argon,  krypton,  and  xentm.  The  stated  fluence  requirements  in  each 
I^boton  energy  tun  was  converted  to  a  required  source  yield,  first  by  satisfying  die  fluence  target 
uniformify  requirement  (Ihble  Q,  and  dien  using  an  estimate  for  die  ddxis  mitigation  distance 
(Tbble  n).  Comparison  widi  projected  yields  fomi  die  JUMTER  design  model  were  given  in 
IhUe  nL  Antidpating  imptovements  in  load  physics  understanding  and  control,  hi^  inqdosiQn 
—  0)  ~  Via]  as  well  as  low  implosion  (~  Vs)  cases  were  reported.  In  die  lower 
energy  ttins,  A  (0- 1  keV)  and  B  (1  >5  keV),  the  calculated  yidds  firomeidier  conqiression  case 
are  close  to  or  exceed  those  required.  This  reflects  die  fact  diat  for  argon,  JUPTlERinqilosions  can 
reach  tte  efficient  seating  reginie  where  die  yield  scales  with  had  kinetic  energy.  Butindiehi^iest 
bins,  D(15  -  30  keV)  and  E  (30  -  60  IceV),  die  {Resent  calculations  indicate  diat  predicted  yields 
fall  short  of  die  required  ones.  Here,  the  inqilosums  are  in  die  low  17  r^inie  and  die  subsequoit 
yields  are  suscqitible  to  the  details  of  the  stagnation  pl^sics,  such  as  temperature  equilibration  as 
shown  above. 

Both  the  inevious  refiort  and  the  present  one  predict  that  low  Z  elements  will  be  efficient 
radiators  on  a  JUPFIER-dass  generator  intheO-landl-5  keV  {dioton  oieigy  bins.  The 
optimal  K-shell  yields  of  aluminum  (>1.6  keV)  and  argon  (>  3.1  keV)  are  consistently  projected 
with  the  J-scaling  law,  as  well  as  2*IMP  simulations,  to  be  '^'30%  of  the  load  kinetic  energy.  But 
these  reports  caimot  be  viewed  as  die  end  of  dieoretical  research  on  JUPITER  PRS  loads.  Firm 
predictions  for  JUPTIER’s  yield  cipability  above  10  keV  are  not  in  hand  due  to  a  number  of  load 
physics  and  related  simulation  modeling  issues.  Hve  of  these  issues  are  discussed  next 

(i)  The  {Resent  r^iort  noted  that  the  pinch  could  not  be  treated  as  a  plasma  shell  inploding 
through  a  vacuum  onto  the  axis.  It  was  found  necessary  to  include  a  core  plasma  which  heats  up 
and  contributes  to  the  softening  of  the  primary  load’s  collapse  on  axis.  Exfierimental  evidence  for 
early  pifl^^a  formatirni  on  axis  was  noted  by  Aivazov  [Ref.9]  for  aluminum  wire  loads.  These 
authors  also  suggested  that  the  blowoff  from  each  wire  during  the  initial  current  rise  subsequendy 
collipses  prior  to  the  main  plasma  mass  due  to  its  low  density.  Other  theoretical  analyses  by 
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Bloomberg,  et  al.,  [Ref.25]  led  to  the  conclusion  that  one  does  not  expect  the  later  evolution  of  the 
wire  plasma  to  be  dependent  upon  initial  conditions.  However,  a  re-investigation  of  this  problem  by 
Whitn^,  et  al.  [Ref.26]  using  a  similar  approach,  found  that  SATURN  kilovolt  yields  vary  directly 
with  die  theoretical  ten^ieratures  achieved  in  the  early  wire  eqdosion  phase  of  the  pinch  dynamics. 
Indirect  evidence  also  exists  for  the  inqportanoe  of  a  pre-fonned  core  plasma  and  enhanced  tran^Kitt 
coefficients  for  argon  gas  puffs  [Ref.lO].  The  sifflulatioos  in  this  lefereooe,  v^iidi  were  matched  to 
esq)eriiiiratal  pinch  densities  and  tenqierahiies  at  implositm,  surest  diat  die  primary  load  bounces 
ofifahotcoceradierdiancdhqisingallthewi^todieaxis.  The  origin  of  die  oofefdasma  for  gas 
puffs  r^ects  the  inidal  density  distribution  arising  fitom  die  nozzle  jet.  Le„  die  r^gioo  within  die 
gaspuffsheUisnotavacomiL  The  wire  arr^  diagnostics,  the  detailed  analysis  of  argon  gas  puffs, 
and  die  present  stu^  to  match  a  large  set  of  SATURN  data  all  point  to  the  potendalty  important 
irifloeaoe  of  die  startup  phase  for  pinch  (fyriainics.  Theteq»oseoftiieearfywireUowofforgas 
puff  density  inofile  to  die  larger  dl/ift  of  JUNTER  compared  to  present  twrfiitw^  is  presentty 
unknowiL  These  aqiects  need  detailed  inwesdgadon.  eqiecialty  rince  a  pre-formed,  axial  plasma 
would  hinder  the  potential  forhi^cmrqxession  [.Ro(<  s  Oj/JRsCirmy)  ~  10]  pindies.  Based  on 
die  LIA  design  in  diis  rqiort,  sudi  pindies  were  found  to  be  essential  if  one  hopes  to  matA  the 
radiadcm  requirenoents  in  the  photon  energy  Uris  above  S  keV. 

(ii)^ln  order  to  achieve  coinous  amounts  (>  Ml)  of  bigb  energy  i^iotons  (>  10  keV)  it  will  be 
necessary  to  inqilodehi^  atomic  number  (Zn«c  >  18)  material  to  suffilden^  high  vdodty  that  it 
is  ioiuzed  into  die  K-shell  stage  upon  diennalizadmi.  For  JUPTIER  designs  whidi  reach 'vdOMA 
peak  load  currmt  this  means  diat  the  irrqilosions  will  not  have  oiough  Idnedc  energy  to  themtalize 
the  bulk  of  a  high  Z  load  into  the  K-shell.  or  it  will  not  cany  sufficient  mass  to  radiate  effidendy.  In 
the  parlance  of  RefJ.c^itinaal  radiating  loads  cm  JUPITER  have  17  'w  1  for  kryptmi.  and  17  <  1  for 
xenon.  The  presmit  ZPIMP  study  of  krypton  <»  a  JUPTlER-dass  generate  clearly  demonstrated 
the  sensitivity  of  krypton  K-sbell  yidds  to  modoate  changes  in  the  mean  plasma  density  and 
tenqierature  at  iiiq>losi<Hi.  It  should  also  be  obvious  that  the  high  energy  radiadmi  yield  will  be 
susceptible  to  plasma  gradients  at  inqilosion,  not  only  because  distinct  r^ons  of  the  plasma  will 
radiate  differmtiy  than  the  average  conditimis,  but  also  because  radiatirm  transput  dqiends  iqxm 
local  gradients.  The  potential  for  stability  and  enhanced  yields  dirough  innovative  load  designs 
can  be  ascertained  cmly  if  die  details  of  the  plasma  dynamics  are  properly  treated.  This  includes 
a  better  understanding  through  ejqferiments  and  modeling  of  the  instabilities,  zippering,  enhanced 
transport  coefficients,  and  general  turbulence  acting  in  the  plasma  to  soften  the  implosion  [Ref.S] 
and  broaten  the  radiation  pulse  width.  The  present  phenomenological  treatment  of  this  stagnation 
{dtysics  through  the  viscosity  paiameto'  tied  to  SATURN  data  makes  the  extension  of  the 
model  to  JUPITER  load  currents  uncertain.  Research  is  under  way  to  see  if  the  comparisons 
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with  SATURN  made  in  Fig.4  lead  to  similar  viscosity  models  for  the  other  existing  pulsed  power 
generators:  PHOENIX,  BLACKJACK  S,  DOUBLE  EAGLE.  Altogether,  these  comments  nialfi» 
the  case  for  load  studies  with  multi-zone  1-D  and  2-D  simulations  which  concentrate  on  the  plasma 
prafites  during  iix4>losion.  Initial  wtxk  on  2-D  stability  has  begun  [Ref.27],  but  it  is  a  large 
undertaldiig  for  die  details  dqieod  cm  the  histxxy  of  die  pinch  when  coupled  to  the  driver. 

(iii)  Not  only  is  die  plasma  dynamics  in^xxtant  in  hi^  Z,  low  17  implosions  for  the  hi^  energy 
photon  yields  meittitxied  in  (i)  and  (ii)  above,  so  is  die  momic  physics  and  ionizption/ejccitatum 
dynamics.  For  hig}i2^«e  radiators,  the  pre-in^dosioa  radiative  losses  through  L- and  even  M-diell 
emissioa  10b  fipom  die  iiqdosion  energy  and  tfaer^  lower  die  iofiizadon  level  at  thermalizadon. 

A  ayrtematie  iiiweWigatinii  iiiln  the  byptftn  jwrid  dependenee  on  the  mnHiJing  nf  the  I  ..Aril 

presented  by  Davis  (Ref28]  at  die  second  JDOST  meeting.  Both  the  krypton  DCA  model  used 
in  the  present  work  and  die  one  in  Ref.23  should  be  amended  widi  an  extensive  M-shell  model 
Improvements  to  die  atomic  i^^sics  model  for  krypton  are  strai^itforward,  but  time  consuming 
in  code  iqidating.  For  xenon,  the  potential  pre-K-shell  losses,  as  well  as  die  the  atmnic  pl^sics 
details,  are  more  severe.  B  has  been  custmuoy  to  me  an  average  atom  (AA)  to  treat  very  high 
Z  material.  It  offers  the  ad|vantages  of  rqiid  calculation  ahile  at  the  same  time  smne  treatment 
ofoorrqdexl^  andM-sheUelectiODicstinctoie.  Hor  rever,  because  various  rates  are  scaled  fimn 
hydrogenic  formulas  it  also  tni^  overeoqihasize  the  effect  of  radiative  transidons.  Thus  the  poor 
yield  from  xenon  in  IhUe  m  nu^  be  indicadve  of  die  difficulQr  in  adhieving  die  K-shdl  ionization 
stage  due  to  die  dominant  Lrand  M-shell  energy  sinks,  or  it  mqr  be  cmly  a  leflecdtm  of  the  inhemit 
assunqiidons  of  die  AA  qiproadL 

(iv)  Closely  omnected  to  die  itnizadcm  (fynamics  of  hi^  Z  materials  is  the  radiadtni  transport 
of  the  free-bound  continuum  emissioiL  IntheJ-scalinglawofReLS,  dietotalK-shellemissimiis 
included  in  yield  esdmates.  For  the  present  wmk.  ZFIMP  uses  a  probdrility-of-escape  tedinique 
f(H' the  fine-bound  emissitm  which  transports  all  the  emitted  radiation  as  if  it  had  the  frequency  at 
die  edge.  In  mder  to  increase  the  yield  in  the  highest  energy  (>  30  keV)  bins,  load  conqxisition 
and  design  may  be  engineered  to  put  out  most  of  the  high  energy  radiadtm  from  recombination 
to  K-shell  states.  In  diis  case  a  nudti-frequency  radiation  transport  is  die  tedmique  for  properly 
calculating  the  emitted  qiectrum.  Since  the  He-  and  H-Iike  continuum  for  kzyptmi  and  xenoi; 
extmds  far  b^ood  the  edge,  the  opacity  evaluated  at  the  edge  is  inaccurate.  In  tlw  high  energy 
jriioton  Inns  the  shqie  of  the  spectrum  is  irxqiortant  for  determining  what  is  essentially  a  small 
component  of  the  total  emitted  radiation.  One  should  also  note  that  by  the  nature  of  an  AA  atomic 
(riiysics  model,  the  free-bound  continuum  from  the  K-shell  stages  is  underestimated  whenever  the 
mean  charge  state  <  Z  >  is  less  than  Zn«c  —  2.  Because  the  bound  energy  levels  in  a  AA  model 
ate  not  as  n^adve  widi  a  higher  <  Z  >,  the  free-bound  continuum  is  shifted  to  lower  energies 
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than  it  should  be  for  17  ~  1  implosions  whetein  the  bulk  of  the  plasma  is  not  in  the  K-shell.  Thus 
a  realistic  estimate  of  the  highest  eneigy  emission  is  best  acconq>lished  with  a  DCA  atomic  model 
[point  (iii)  above]  in  conjuncti<»  with  multi-frequency  radiation  tranqxMt  The  drastic  conclusion 
Aat  a  PRS  cannot  produce  the  required  high  energy  (>  30  keV)  radiation  yield  and  tme  nuist  look 
to  other  techniques^  is  at  die  least  sinqily  premature.  Future  work  on  the  high  energy  conqxment 
could  address  the  problem  of  xenon  ionization  ttynamics  widi  a  modified  DCA  atimik  model  and 
multi-firequency  transport 

(v)  Rnally,  die  prosentZPIMPoode  is  designed  to  tun  off  of  a  lunqiedciicuit  model  for  die 
dicttit  This  is  acoqitidile  for  the  lineac  voltage  adder  fix  whidi  the  vidtage  pulse  origiiutfes 
away  in  time  ffooLdie  load  and  die  load  behavior  does  not  alter  the  pulse,  excqic  for  long  M  200 
nsec  inqdosion  times.  But  for  any  indnctiue  storage  design,  the  driving  open  dronit  wdtage  is 
effectively  diat  created  across  die  plasma  opening  switch.  Hus  efement  is  dose  to  the  load  and 
die  voltage  does  requmd  to  die  load's  bduwioc.  The  scaling  law  calculations  of  Ref  J  do  treat 
die  ciicnit  as  a  transmission  line,  but  ZFIMP,  and  mote  advanced  multi-zone  1-D  and  2-D  codes, 
should  be  revised  to  do  likewiae.  TUs  would  allow  detailed  nnmlatioomodds  to  address  qiedfic 
aspects  of  power  flow  and  load  ooupliiig  on  varkms  madiine  designs,  particolatty  if  die  revised 
code  also  included  electron  losses  in  magnetically  insolatnd  transmission  lines.  The  cadra  difficulty 
in  setf-consisteod^  coupling  die  eneigy  fiom  a  transmission  line  circuit  into  die  front  end  load 
region  lies  in  die  (fiffusion  of  die  cxteriud  magnetic  field  inside  the  (dasmar^ion.  TbispioUem 
is  presently  under  developmenL  The  revised  STMP  code  should  also  investigate  radiation  yields 
firom  the  recent  JUPITER  deagn  verrions  for  die  llA  and  lES.  Both  deliver  more  kinetic  eneigy 
than  the  circuit  used  in  die  present  pqier  and  offer  die  potential  for  larger  yields  than  LIA-1  in  die 
energy  bins  of  Thble  m.  fri  diese  new  designs  Ref3  showed  that  improved  K-shell  yields  fix’ argon 
could  be  adiieved  if  die  pindi  lengdi  was  increased  form  4  to  6  cm,  vdiile  shorter  pindier  lengths 
(~2  cm)  {Moved  benefidal  for  xenon. 
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